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Abstract 
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approach to computation of the optimal perfomance index achievable by causal pre- 
coders for given channels. This quantity is important in the numerical determination 
of optimal precoders. 
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1 Introduction 



Many of todays state-of-the-art wireless systems adopt multiple-input 
multiple-output (MIMO) transmission to increase spectral efficiency together 
with multi-carrier methods to cope with intersymbol interference (ISI). Multi- 
carrier methods simplify channel equalization because they decompose fre- 
quency selective channels into multiple flat fading channels (the so called car- 
riers), which can be easily equalized. While multi-carrier transmission offers 
many advantages including effective channel equalization, it also exhibits some 
drawbacks regarding the peak-to-average power ratio (PAPR). Often single- 
carrier transmission, where the frequency selective channel is approached di- 
rectly, is considered as an alternative to multi-carrier transmission [Tf2~f3] . 
While therefore single-carrier transmission is interesting on its own, it has been 
further shown in [4f5] . that in fact many common multi-carrier, code-multiplex 
and space-time block-code systems can be modeled as single-carrier systems 
by virtual enhancement of the MIMO system. Various authors used this ap- 
proach to derive new equalization methods based on single-carrier equalization 
in order to exploit joint equalization of spatial, time and code or frequency 
domains [4f5f6l[7] . There, and generally for linear time-invariant (LTI) equal- 
ization of single-carrier systems with zero- forcing and causality constraint, one 
usually solves the so-called Bezout Identity 

H(e ie )G(e ie ) = 1 (0 < 9 < 2vr), 

where the matrix- valued transfer function if of a stable and causal LTI system 
(the frequency selective MIMO channel) is given, and a transfer function G 
of a stable and causal LTI precoder, which equalizes H, has to be computed 
[8]. Transmitters may use such G to pre-equalize the channel. Alternatively, 
receivers can also solve the Bezout Identity for the transposed channel H T (i.e. 
H T G = I) and equalize the channel H with the transposed solution G T . The 
main difficulty in solving the Bezout Identity is the causality of G, because 
the naive approach 

G(e w ) = H(e ie y{H(e ie )H(e ie Y}- 1 (0 < 6 < 2vr) 

of a pseudoinverse generally results in a non-causal precoder [9]. If the number 
of channel inputs equals the number of channel outputs, the pseudoinverse is 
the unique solution to the Bezout Identity. The situation changes if the num- 
ber of inputs of H is larger than the number of outputs. Now precoders for H 
no longer have to be unique. Usually this non-uniqueness then is exploited to 
choose a causal G that is optimal in some sense. The two common optimal- 
ity conditions are minimality of the equalizers energy and minimality of the 
equalizers peak value, respectively. The minimal energy condition corresponds 
to the classical approach of signal-to-noise ratio (SNR) maximization [5f7lfT0] . 
However, this approach is only feasible if the statistical properties of the noise 
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are known. For unknown noise statistics, it cannot be applied. Picking up an 
idea from robust control (see e.g. [11]), where one is concerned with unpre- 
dictable errors that arise e.g. due to uncertain modeling, instead minimization 
of the equalizers peak value has been proposed [6J. As we shall see later, this 
minimizes the worst case error instead of the average error, which cannot be 
determined due to the unknown noise statistics. 

In this paper we are interested in the optimal performance that causal pre- 
coders and equalizers can archive regarding the worst case error. Therefore 
we show how a solution to the Bezout Identity with minimal peak value can 
be constructed. We discuss why this gives the best upper bounds on various 
perturbations in the system. Contrary to other ways to solve the Bezout Iden- 
tity, our construction holds for the most general case of systems with infinite 
impulse responses (which are not required to be rational) and even infinite 
input and output vectors, i.e. we allow systems to have infinite temporal as 
well as infinite spatial dimension. We further give a new result on the numer- 
ical computation of the minimum peak value achievable by causal solutions 
to the Bezout Identity if the numbers of inputs and outputs are finite. This 
is important because for all methods known to the authors that solve the Be- 
zout Identity with minimal peak value in a numerically exploitable way, the 
minimal peak value has to be known in advance [6|12| . Therefore efficient com- 
putation of the minimal peak value is important for numerical solution of the 
Bezout Identity. We point out that the optimization approach in [13j requires 
no prior knowledge of the minimal peak value. However, it only computes 
finite impulse response solutions to the Bezout Identity, which are generally 
suboptimal. 

We proceed as follows. In Section 2 we give our problem statement after we 
have introduced some notation and necessary basic mathematical concepts. 
We further discuss the practical interpretation of our problem statement. In 
Section 3 we derive our results on the numerical computation of the minimal 
peak value achievable by causal solutions to the Bezout Identity. Then a op- 
timal causal precoder is constructed in Section 4. We finally draw conclusions 
in Section 5. 



2 Preliminaries 

2.1 Notation 

We denote the complex numbers by C, the complex matrices with m rows and 
n columns by C mxn and the complex column vectors by C m := C mxl . The 
complex unit disc is given as D := {z E C: \z\ < 1}, its border is the unit circle 
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T := {z G C : \z\ = 1}. Complex conjugation is denoted by (•), taking adjoints 
in a Hilbert space by (■)*. Furthermore H,£ and £* denote separable Hilbert 
spaces with scalar products (■, -)^, (-, -)e and (■, -)^, respectively. By Ti® £ we 
mean the direct Hilbert sum, i.e. the space TixS equipped with scalar product 
(h(B e,g f)n®£ := (K g)n + ( e , f)e- The space of bounded linear operators 
between £ and £* is denoted by £(£,£*). It is equipped with the operator 
norm ||T|| op := sup e6 £ i| e || £: =i HFe]^. On any space the identity operator is 
written as I. For matrices A G C mxn the smallest and largest singular value 
will be denoted by a m i n (A) and cr max (A), respectively. The closure of a set M 
is denoted by closure M, the space spanned by all linear combinations of its 
elements by spanM. 



As usual, Lj(X) denotes the space of (equivalence classes of) p-integrable 
functions on T with values in a Banach space X. The norm in Lj(X) is 
11/11? : = Ie=o\\f(^)\\x^ for 1 < V < oo and \\f\U := esssup feT ||/(C)||x = 
inf {m > : yu({C £ T : ||/(C)l|x > m }) = 0} for p = oo, where // denotes the 
Lebesgue measure. We refer to [HI Section 3.11] and the references therein 
for details on integration of vector- and operator-valued functions. If p = 2, 
Lj(£) equipped with the scalar product (/, g) 2 ■= Jelo{f(e ie ), g(e l6 )) £ ^ is a 
Hilbert space. For F G we denote the point-wise adjoint by F*, 

i.e. F*(() = (F(())* almost everywhere on the unit circle. 



2.2 Basic Results and Concepts 

2.2.1 Hardy Spaces and Toeplitz Operators 

We introduce the usual Hardy spaces on the disc by 

f r 27T d9 1 

Hl{£) := lu : D -> £ : u analytic, \\u\\ 2 2 := sup / \\u(re ld ) < oo }, 
[ o<r<iJe=o 2it J 

H™(£, £») := { F : D -> £(£, £,) : F analytic, := sup || op < oo 1 . 

t zeo J 

The Hardy spaces play an important role in systems theory, since they are 
the set of transfer functions of causal finite energy signals and causal and 
energy-stable transfer functions for LTI systems, respectively pj)]. Definition 
is also possible on the upper half plane instead of the unit disc. On both 
domains, the Hardy functions are completely determined by their values on 
the borders of the domain. Therefore, each Hardy function on the unit disc 
has a corresponding function on the circle. The space of those corresponding 
functions can be given as follows. 
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For functions / G L\{£) or / G L^(£(£, £*)) the fc-th Fourier coefficient is 

/,:= r /(e^)e-^f (* G Z). 
Therewith, the Hardy spaces on the circle are given by 

:= {u G L£(£) : M fc = for k < o} , 
H?(E,£*) := {F G L$(C(£,8.)) : F fc = for fc < o} . 

It is important to know that the two notions of Hardy spaces on disc and 
circle are equivalent, since by Fatou's Theorem the radial limit (bu)(e ld ) : = 
lim r yi u(re 10 ) exists almost everywhere and the mapping b is an isometric 
isometry between the Hardy spaces on disc and circle (see [331 Th. 3.11.7, 
3.11.10]). Therefore we will only explicitly distinguish between those spaces if 
necessary, and simply write H 2 (£) and H°° (£,£*) otherwise. 

An important property of L\{£ ) is ParsevaVs Relation ([HI p. 184]), by which 

oo 

ll M ll2 = for all it G L|(£). 

k=— oo 



We will now introduce Toeplitz operators, which are the standard example for 
operators on Hardy spaces and which will also play an important role in what 
follows. The orthogonal projection (P + u)(Q := Y^k=o^kC k from L\{£) into 
Hj(£) is called the Riesz Projection. The projection from H%(£) into the space 
of degree N polynomials is (P N u)(C) := Ek=o^k( k - Now for F G Lf(£,Q 
the Toeplitz operator with symbol F is the operator which maps Hj(£) into 
H%(E*) via Tpu : = P+(Fu). 

The next result allows us to get an exact estimate of the minimum norm 
achievable by solutions of the Bezout Identity. 

Theorem 1 ([H, Th. 9.2.1]) Let F G H°°{£,£*) and 5 > 0. Then some 
G G H°°(£*,£) with \\G\loo < and F{z)G{z) = I for all z G © exists if 
and only if 

\\T F *u\\2 > S\\u\\ 2 for all u G H 2 (£*). 

2.2.2 Schur Class 

Functions in the unit ball of H°°(£,£*), the so-called Schur class 
S (£,£*) :={F E H°° (£,£*) : H^Hoo < 1} , 
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have some special properties, which will turn out to be useful in the construc- 
tion of a minimum norm right inverse. Every Schur function can be factorized 
as follows. 

Theorem 2 ([IH Th. 2.1]) Let F : D -> £(£,£*). Then F e S{£,Q if 
and only if there exists a holomorphic function W : © — > £(7Y,£*) such that 

I - F(z)F(w)* = (1 - zw)W(z)W(w)* (z, w G D). 

Note that W can be given explicitly, see [TTJ, Sec. 3.3]. We finish with the 
observation that also certain block operators define Schur functions. 

Lemma 3 ( [12l Lem. 2]) Let T e £(H®£,H®£*) with \\T\\ op < 1. Tften 
T /«as a unique block matrix representation 



A B 




n 




n 


C D 




£ 




£ 



and the function 

F : O — >• £(£,£*), F(z) := D + Cz(I - zA)- x B 
is Schur, i.e. F <G S(£,£*). 

Functions defined as F in the Lemma above are known in operator theory as 
characteristic functions, while unitary operators like T are known as unitary 
colligations. Those concepts resemble much the concept of a transfer function 
and a state-space realization in control theory. We refer to |T7l[T8] for details. 

2.3 Problem Formulation 

Before we give an exact problem formulation we introduce and discuss the 
target objective 

lopt (H) := inf ({HGIU : G G H°°(£., S), H(z)G(z) = J for all zeB}U {oo}) , 

which is, as we will see, a tight lower bound on the worst-case transmit energy 
enhancement of causal precoders for the channel H, and a measure for the 
achievable robustness against imperfectly known channel transfer functions. 
Note that in particular 7 opt (i7) = oo if and only if H has no right inverse in 
H°°. We always assume H e H°°(£, £*) unless we explicitly mention otherwise. 

It was shown in [19] that if dimf* < oo, existence of a right inverse in H°° is 
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further equivalent to 

H(z)H(z)* > 5 2 I for some 5 > and all zeB. 

It is somewhat surprising that although by the result from [19J 7 op t < oo if 
and only if 

8 C := sup [5 > I H(z)H(z)* > 5 2 I for all z £ d| > 0, 

<5 C has no direct connection to 7 op t, i.e. 7 op t cannot be computed from 5 C 
[9]. However, as we will see, it is important to know 7 opt in advance of the 
construction of an optimal precoder. Therefore we derive a new method for 
numerical computation of 7 opt and then solve the following problem. 

Problem 4 Let j op t{H) < 00. How can G E H°°(£„£) with H{z)G{z) = I 
for all z € IB) and ||G||oo = Japt(H) be constructed? 

We close this section with a short discussion in which sense minimization of 
the infinity norm in Problem 0] gives optimal filters. The input-output relation 
of a frequency selective MIMO channel is given by 

y(0 = H(()x(0 + n(0 (C£T), 

where H denotes the channel, x the transmitted signals and y and n the 
received signals and additive noise, respectively. If a precoder G for H is used 
to pre-distort the transmitted signals, this input-output relation changes to 

y(C) = tf(CMCMC) + n(C) = *(0 + n(C) (C e T). 

There are two advantages in minimizing the infinity norm of the precoder. 

The first advantage is minimization of the transmit signals energy. The energy 
necessary to transmit a signal x using the precoder G is given by 1 1 C?a? 1 1| . 
Without loss of generality, let us assume that = 1. Then, it can be 

shown that the transmit energy necessary in the worst case is exactly || C|| ^ , 
i.e. 

sup ||Cx|| 2 2 =||G||^. 

x&H 2 {£ t ),\\x\\ 2 =l 

Thus, minimizing ||G||oo guarantees the lowest amount of necessary transmit 
energy. If equalizers instead of precoders are considered, i.e. 

y(0 = G(C)[tf (CMC) + "(C)] = x(C) + G(C)n(C) (C G T), 
this is equivalent to minimal worst case noise enhancement. 

The second advantage of minimization of the infinity norm is robustness. As- 
sume an imperfectly known channel transfer function if a = H + A with right 
inverse Ga, where H is the correct channel and A is a perturbation. Using 
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the same argument as before, we see that the energy of the worst error that 
can result from the perturbation equals 

sup \\x — HGax\\1 = sup || ACa^III = I|AC?a||to- 

x£H 2 {£»),\\x\\ 2 =l xeH 2 (£,),\\x\\ 2 =l 

Since it holds ||AGa||^o < II AH^HGaII^ and this inequality can become sharp 
e.g. for A = 51, we see that minimizing ||Ga||oo also minimizes the worst case 
error that results from an imperfectly known channel transfer function. This 
argument applies to equalizers in the same way it applies to precoders. 



3 Computation of the Optimal Norm 

This section deals with the computation of the optimal norm 7 opt achievable 
by solutions to the Bezout Identity. Since many algorithms which directly 
solve Problem [4] only compute suboptimal solutions, i.e. given 7 > 7 opt they 
compute a right inverse G 7 with norm ||G 7 || 00 < 7, it is important to know the 
optimal value for 7 in advance [6lll2| . We point out that computation of 7 opt 
also arises in other contexts, see e.g. Remark 1 in [20] (with the next corollary 
in mind). 

We start with an exact (but incomputable) formula for 7 opt . The next two 
corollaries are direct consequences of Theorem CD 

Corollary 5 For p(H) := inf w6 H 2 (f»).ll«ll2=i \\ t hM\2, it holds j op t(H) = 
p(H)-\ 

Corollary 6 If j opt (H) < 00, a right inverse G E H°°(S^,S) with ||G||oo = 
7 opt (#) exists. 

The interesting thing about Corollary [5] is that it shows us why the optimal 
causal equalizer cannot perform better than the optimal non-causal one. Note 
that the optimal norm for non-causal equalizers is given by 

inf \\H*u\\ 2 \ 

ue// 2 (£*),|M| 2 =i / 

(see [9]), which is the same formula as Corollary [5l except for the additional 
Riesz projection P + : 

7**iH)=( mf ||P + (^)|| 2 ) . 

It is now clear that causal equalizers perform worse because the signal energy 
of u which is mapped into the non-causal part of H*u is cut off. How much 
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energy is shifted into the non-causal part thereby depends on the Fourier 
coefficients of H*, which are related to H by H* k = H*_ k for k G Z. 



We now derive a computable approximation of 7 op t. The main idea will be to 
approximate the relation 7 opt = p^ 1 from Corollary El In order to compute 
7o P t, we try to approximate p with 

p N (H):= inf \\P N T H *u\\ 2 , 

«eP]V-H^(t*),||«||2 = l 

i.e. we restrict domain and image of T H * to polynomials of degree N and 
take the infimum for this restriction. Because PnTh*Pn is linear and finite 
dimensional, it can be represented by a matrix. 

The main result of this section is the following. 

Theorem 7 The sequence {pn(H)}n is monotonically decreasing and con- 
verges with limit 

lim p N {H) = p{H) = lopt {H)- 1 . 

N— >oo 

IfH e H°°(C mxn ) with m<nQ and 



r 



H,N 



H*H{... H 



N 



H* 



H 



N-l 



... H* 



£ ^n(N+l)xm(N+l) 



Pn can be computed as Pn(H) = o- min (TH^). 



PROOF. We only sketch the proof here, the full proof is given in the ap- 
pendix. It consists of three main steps. The first step is to show that the 
sequence {pn(H)}n is monotonically decreasing and lower bounded by p(H). 
The main idea is that the relation 

Pn(H) = inf \\PnTh*u\\2 = inf llTViilU 

u£P N H 2 (£«),\\u\\ 2 =l u£P N H 2 (e*),\\u\\ 2 =l 

holds for every N E N and thus the infimum is always taken over the same 
target objective, but over a space which increases with N. This is done in the 
appendix in Proposition [T4l In a second step it is shown that the lower bound 
p(H) for {p N (H)} N is sharp. Therefore for arbitrary e > a sequence {u^}n 
such that 

u N e PnH 2 (£*), \\un\U — 1 an d hm \\PnTh*u n \\ 2 < p(H) + e 

N— too 



1 Note that trivially ^y opt (H) = oo for m > n. 
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is constructed in the appendix in Proposition [TH1 Thus Pn(H) converges to 
p(H), which is equal to ■y opt (H)~ 1 by Corollary [5l Finally Proposition [TBI in 
the appendix gives the formula for computation of Pn(H) via singular value 
decomposition if H is matrix-valued. ■ 

Since the arguments used to prove Theorem [7] hold analogously if we approx- 
imate 

sup ||Tff»w|| 2 = ||Tff*|| op = ||T^|| op = ||Tf/|| op = ||-H"||oo 

u£H 2 (C m ),\\u\\ 2 =l 

instead of Pn(H) = mf u€H 2( Cm ) ,|| u || 2 =i ||Th*w||2, we also see that for H G 
i7°°(C mxn ) the sequence {c mauX (TH,N)}N is monotonically increasing and con- 
verges with limit 

lim cx max (r^Ar) = \\H QQ. 

We note that the well-known fact that the limit ||-H"||oo of a max (T HN ) can be 
found by performing a grid search over all frequencies, i.e. 

lim cr max (T H ,N) = \\H\\oo = esssup CeT o- max (i?(C)), 

does not carry over to computation of 'y op t(H). Here, in general we have 

lim 0- m i n (r ff) jv) = 7o P t(^) _1 7^ essinf CeT a min (if(C)). 

N—>oo 

This dichotomy results from the fact that while indeed 

sup ||f/"*it|| 2 = sup \\P + (H*u)\\ 2 , 

u£H 2 (S t ),\\u\\ 2 =l u£H 2 (S,),\\u\\ 2 =l 

in general we have 

inf \\H*u\\ 2 ^ inf \\P + {H*u)\\ 2 . 
This can be easily seen in the next example. 

Example 8 Set H(() = ( for ( G T. Then by Parseval's Relation 

inf ||i/*-u|| 2 = inf ll-ulU = 1, 
«e.ff 2 (c),|M|2=i «G// 2 (c),||n|| 2 =i 

however for u(z) = 1 we have (H*u)(() = ( and therefore 

\\P + (H*u)\\ 2 =\\0\\ 2 = 0. 

It is also important to note that Theorem [7] does not generalize to the case 
H G . We give an example where p{H) = 1, a inverse in H°° exists, but 
the smallest singular values of the finite sections do not converge to p(H). 
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Example 9 Set H(() := ( for ( G T. Then by ParsevaVs Relation 



inf ||Th-*m|| 2 = inf II T>w II 2 = inf llnlU = 1. 
ue// 2 (c),||«|| 2 =i «eH 2 (C),|M| 2 =i mg-H" 2 (c),||u|| 2 =i 



Further, H has a inverse in H°°, i.e. G(() = (. However, 





/ 


" H* 


H{ ... 


TJ* 






/ 








\ 


"~min 




Ht x 


H* •-. 


H* 








l '• 








V 


TJ* 

_ —N 


... Hl x 


H*_ 


) 




V 




1 


/ 



= 



for all N eN. 



4 Construction of the Optimal Causal Precoder 

In this section we construct a minimum norm solution to the Bezout Identity, 
i.e. we solve Problem 01 The major idea of the proof is the following. We 
first show how to construct right inverses with norm at most one. Then given 
any H G we apply this technique to the scaled function joptH. 

Appropriate rescaling of the obtained inverse will result in a minimum norm 
right inverse. 

4-1 Schur Right Inverse 

The first step is construction of a Schur right inverse. Therefore we factorize 
the function to be inverted similar to Theorem [2] and use this factorization 
to construct a contraction of the form of T in Lemma El The characteristic 
function of this contraction then is the wanted right inverse. This is a variant of 
the technique known as "lurking isometry method", which has been introduced 
by Ball and Trent [171 Th. 5.2] and independently Agler and McCarthy [21] 
to solve the Bezout Identity. 

We start with the factorization. 

Lemma 10 Let H have a right inverse G G S (£*,£). Then there exits a 
holomorphic function W : D — > C{7i,£*) such that 

H(z)H(w)* —/ = (! — zw)W(z)W(w)* (z,«ieD). (1) 
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PROOF. By Theorem [2] there exists a holomorphic function W : D 
5) such that I - G(z)G(w)* = (1 - Thus 

H(z)H(w)* - H(z)G(z)G(w)*H(w)* = (1 - zw)if (2)W(z)W(w)rff (w)*. 

Since #G = / we obtain with W(z) := H(z)W(z) that 

H(z)H(w)* - I = (1 - 



We can now introduce the appropriate block operator. 

Definition 11 Let if have a decomposition like (DO) in Lemma\T^ We define 
the sets 



D : = closure 



Ro := closure 
and a function Vq : D — > i? % 







span < 




V 




/ 




span < 




V 





u>PF(io)* 
H(wY 



W(w) 
I 



fc=0 



H(w)* 



k=0 



W(w) 

I 



Note that it can be easily shown with (pQ) that Vq is a isometry, i.e. 



Vn 



h 


,v 


h 


> -< 


h 


5 


h 


e 




e 




e 




e 



for all 



eH®£. 



Later we will use this fact when we apply Lemma [3] to an extension of Vq. The 
wanted right inverse can now be given explicitly. 

Theorem 12 Let H have a decomposition like (Dp) in Lemma\T(^and construct 
Vq as in Definition [771 Denote by 



Vr 



00 



A B 




H 




H 


C D 




£ 




£ 
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Vood 



the continuation of Vq with zero, i. e. 

, d i D ' 
Then the function 

G(z) := D* + B*{I — zA*)~ 1 zC* (z G D) 
is a Schur right inverse of H , i.e. G G S (£*,£) and HG = I. 



PROOF. Let w G 


D. By construction of Vqo 


it holds 




A B 




wW(w)* 




W(w)* 




C D 




H(w)* 




I 



for all e* G £*, which is equivalent to 

AwW(w)* + BH(w)* = W(w)* (2) 

and 

CwW(w)* + DH{w)* = I. (3) 

Since 1 1 Vbo 1 1 op < 1 1 1 1 op = 1 because Vq is an isometry, we have ||A|| op < 1 and 
thus ||Au)|| p < 1. Thus I — Aw is invertible, and ([2]) yields 

W(w)* = (I - AH))- 1 BH(w)*. 

Plugging this representation of W{w)* into (111) results in 

Cw(I - Aw)- x BH{w)* + DH(w)* = I. 

Taking adjoints and replacing w by z shows that 



H(z) 



D* + B*(I - zA*Y 1 zC* 



This right inverse is Schur by Lemma [31 

4-2 Minimum Norm Right Inverse 



The extension of Theorem [12] from an upper bound one on right inverses 
to arbitrary bounds is a simple scaling argument. Note that in particular the 
upper bound 7 = 7 opt (if ) is valid due to Corollary [6l and results in a minimum 
norm right inverse of H. 
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Corollary 13 Let j op t(H) < 7 < 00. Denote by G G S(£*,£) the right in- 
verse to H := as given by TheoremUM Then G := 7G is a right inverse 
of H with || C|| oo < 7. 

PROOF. Since 7 op t(#) < 7 < 00, a right inverse G G H°°(£*,£) of H with 
|| G\\ oo < 7 exists by Corollary El Thus 

jH 1 - 1 G = I, h^GIU < 1, 

which shows that H = 'jH has a right inverse in S(£*,£). Let G G S{£*,£) 
denote the right inverse of H given by Theorem [121 Then G = 7G holds 

II^Hoo = tII^IIoo < 7 as well as 

EG = -f- l H-fG = I. 



5 Conclusions 

In this paper we considered the problem of the construction of a causal pre- 
coder with optimal robustness for a stable and causal LTI system with multi- 
ple inputs and outputs. This problem is equivalent to finding a solution to the 
Bezout Identity with minimized peak value, for which we gave an explicit con- 
struction. We derived a novel method for numerical computation of the lowest 
peak value achievable in this problem, because it has to be known prior to the 
construction of the optimal precoder. This method is based on computation of 
a singular value decomposition of the finite section of a certain infinite block 
Toeplitz matrix, which is directly constructed from the Fourier coefficients of 
the systems transfer function. 



Appendix 

The complete proof of Theorem [7] follows splitted in three propositions. 

The first proposition shows that {pn(H)}n is monotonically decreasing and 
converges with a limit not lower than p(H). 

Proposition 14 It holds 

Pn(H) > PN+ i(H) > p(H) 

for all N G N. 
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PROOF. Let u G H 2 (£*). We set v := P N u and w := T H *v. A simple 
computation shows that the Fourier coefficients of w = P + (H*v) are given by 



, k < 



Since by construction Vk = for k > N, we see that Wk = for k > N. Thus 



N 



\T H *P N u\\ 



\w\ 



Yl W^k\\l = Yl = \\Pnw\\1 = \\P N T H *P N u\\ 



(4) 



k=0 k=0 

holds by Parseval's Relation for every u G H 2 (£*). 

Because trivially P N H 2 (£*) C P N+ iH 2 (£*) C H 2 (£*), we obtain with (j4j), that 



p N (H)= inf ||Pjv^fr*tt||2 

«ePjv^ 2 (£.).ll«ll2=i 

= inf ||7if.u|| 2 

«eP]vH 2 (^»,),||«||2=i 

> inf ||Tff*d| 2 

ueP N +iH 2 (£*),\\u\\2=l 

> inf IITh.mIU 

ue-H" 2 (£*),||it|| 2 =l 

= p(H). 



Pn+i(H)) 



We now ensure that the limit of {pn{H)}n also is not greater than p(H). 
Proposition 15 For every e > there exists K G N smc/z £/ia£ 

p N (H)<p(H) + e 

for all N> K. 



PROOF. We assume H ^ since the case H = is trivially true. Let e > 
and choose u G H 2 (£*) with ||m|| 2 = 1 such that 



\T H .u\\ 2 -p(H)\ 



\T H *uh ~ „ inf ||Th-.m|| 2 

ueH 2 {£,),\\u\\ 2 =l 



e 

< -. 
~ 6 



(5) 



Since u G i/ 2 (£*), Th*u G H 2 (£) and ||w||2 = 1, Parseval's Relation shows 
that 



lim IIPjvU — mIU = lim IIPaj-Th-.u — T^wlU = 0, lim llPvttlU = 1. 

JV^oo JV-*oo N^oo 
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Thus K G N exists such that 



\\Pnu-u\\ 2 <-\\t h 4; p \ (6) 



\\P N T H *u- T H *u\\ 2 < - and (7) 

6 

\\P N uh> P -^i (8) 



6"~" llop ' 
e 

6 

P 

p(H) 4 f 

for all N > K. 

Then for iV > it follows that 



IIPatTh.Pjvm - T H *u\\ 2 < \\P N T H *(P N u - u)\\ 2 + \\T H *u - P N T H *u\\ 2 



< || iVTjf. || op ||Pjvm-m|| 2 +\\T h *u- P n T h *u\\ 2 



<\\Th* Hop <e/(6||T K , ||op) by {6J 



<e/6 by 



< 



(9) 



and therefore 



\P N T H *P N u\\ - inf ||T H *n|| 2 
«e.ff 2 (£*),|Mj2=i 



<\\\P N T H *PNuh- \\Th**>M 

<e/3 by JU 



||Pff*«||2- ,inf ||Pff*«||2 

«e J H" 2 (£,),||«||2=i 



<e/6 by (© 



< 



We see that 



|PvT„*Pv«|| 2 < inf ||Tfl.ri|| a + - = p(P) + -. (10) 

uEH*(£*),\\u\\2 = l Z Z 



Since || P/vttUa > for N > K by (jHJ), the sequence {mtv}/^^ given by 

Pv« 



l|Piv«lh 



G P N H 2 (£*) 



is well-defined. We obtain the intended result 
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Pn(H) 



- inf 

ueP n H 2 (S*),\\u\\ 2 =: 

< \\PnTh*Un\\2 
\\PnTh*Pnu\\2 

\\Pnu\\2 
(by CSJ) p(H) + § 

- \\P N u\\ 2 

(by ®) 



for all N > K. 



We know now by the Propositions [14] and [15] that the sequence pjv converges to 
p for N — > oo. However it is still unclear, how pjy can be computed explicitly. 
The next proposition gives a simple formula for the numerical computation of 
Pn- 

Proposition 16 Let H G if°°(C m><ri ) with m < n and set 



Hq HI ... H 



N 



H*.. 



H 



N-l 



... H* 



^ £m(N+l)xm(N+l) 



Then p N {H) = a min (T HtN ). 



PROOF. Let USV* = T^ Tv denote a singular value decomposition of r# jv 
with singular values 

0i > • • • > TO (iv+i) > 0- 
Then U G C n ^ N+1 ^ N+1 ^ and V G C m ( iV+1 ) xm ( iV+1 ) are unitary matrices and 

£ g C n(7V+l)xm(7V+l) j g of thg form 

0"! 



0m(AT+l) 



Let u G P/viJ 2 (C n ) and set t> := PnTh*u. We saw already in the proof of 
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Proposition [HI that the non-zero Fourier coefficients of v are uniquely deter- 
mined by the relation 















V N 





H* H* 

o m 



u* 

N 



H 



N-l 



... H* 













Uq 




Uq 




«1 




Ui 




— ^H,N 




UN 




UN 



Thus by Parseval's Relation 



p N (H)= inf \\PnT h *u\\ 2 

«GP]V-f/ 2 (C m ),||M||2 = l 

= inf ||rV iA rw|| 2 

li eC m ( JV + 1 ),||«|| 2 =i 

= inf II Su II 2 

uec m ( JV + 1 ),||u|| 2 =i 

= Cm(AT+l) 

= CTminCr^Ar). 
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